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We report the higher order moments of the net-charge multiplicity dis- 
tributions for the Au+Au collisions at -\/%n = 

7.7, 11.5, 19.6, 27, 39, 62.4 
and 200 GeV in the STAR experiment at the Relativistic Heavy-Ion Col- 
lider (RHIC). The energy and centrality dependence of higher moments and 
their products (such as <J 2 /M, Sa and k<t 2 ) are presented. The data are 
also compared to Poisson expectations and Hadron Resonance Gas model 
calculations. 

PACS numbers: 25.75.-q,12.38.Mh 

1. Introduction 

At extreme high temperature and/or chemical potential, hadronic mat- 
ter undergoes a phase transition where color degrees of freedom plays im- 
portant role, known as Quark Gluon Plasma (QGP). The QCD based model 
calculations claim that at vanishing baryon chemical potential (/xb) the tran- 
sition is simply a crossover [1] . At large //b > the phase transition is of first 
order [21 IH [5] • The presence of these two types of transition, at two extreme 
limit of fj,B, urges to have an end point in the first order phase transition in 
QCD phase diagram, known as the QCD Critical Point (CP). To probe the 
CP and map the phase boundary in phase diagram, the beam energy scan 
program [6l [7] has been undertaken with successful data taking in year 2010 
and 2011. The STAR experiment has collected data for Au+Au collisions 
at ^/i NN = 7.7 to 200 GeV covering large range of /i B = 20 to 410 MeV [8j 
in the phase diagram. 

The various QCD based models reveal, at the CP, the susceptibilities 
and correlation length (£) of the system diverge [5]. The higher order sus- 
ceptibilities and higher power of the correlation length amplify the signal of 
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divergence at the CP. The first four moments, such as mean (M), standard 
deviation (a), skewness (S) and kurtosis (k) of the conserved charged dis- 
tributions like, net-charge, net-proton and net-strangeness are very closely 
related to these respective higher order susceptibilities [5] and higher power 
of the correlation length [9]. For example, the variance, skewness and kur- 
tosis have been shown to be related to the £ 2 , £ 4 ' 5 and £ 7 , respectively [9]. 
In order to cancel the volume term in the susceptibilities, different combi- 
nations of the various moments are constructed, some of these are: a 2 /M, 
Sa and na 2 . These products, as a function of beam energy, are expected to 
show non-monotonic behavior near the CP. Recently in the STAR experi- 
ment, the analysis has been made on the higher moments of the net-charge, 
net-proton and net-kaon [TO] for the beam energy scan program to probe the 
CP in phase diagram. Beside the CP search, recent lattice QCD model cal- 
culations [llj motivate to extract the freeze-out parameter like, temperature 
and chemical potential, from the experimentally measured higher moments 
of the net-charge multiplicity distribution. 



2. Analysis Details 

The STAR experiment has taken data for Au+Au collisions at \/^nn = 
7.7 to 200 GeV in last two years. The STAR detector system provides the 
excellent particle identification and large acceptance for the event-by-event 
fluctuation analysis. The Time Projection Chamber (TPC) is the main 
tracking device. Extensive quality assurance is performed for each energy 
in order to minimize the fluctuation of detector efficiency. 

The event-by-event net-charge distribution is measured for Au +Au col- 
lisions occurring within ±30 cm along the z position of the interaction 
point from the TPC center and 2.0 cm radius in the transverse plane. The 
charged particles are measured between the transverse momentum (pr ) 
range 0.2<p^ <2.0 GeV/c and pseudo-rapidity (rj) range at \r]\ < 0.5 re- 
gion. The standard STAR track quality cuts are used for this analysis. The 
contamination of spallation protons, produced due to beam pipe interac- 
tion, are suppressed by removing protons within 200 <px < 400 MeV/ c. 
To avoid the auto-correlation effect in the higher moments analysis, the cen- 
trality selection has been done by uncorrected charged particles measured 
within 0.5 < |r/| < 1.0 from the TPC detector. The heavy ion collision 
geometry of each event is characterized by the number of participating nu- 
cleons (A^t). The average number of participant ((Ap art )) for each cen- 
trality bin is calculated by the Monte Carlo Glauber simulation. The finite 
centrality bin width effect has been avoided by using centrality bin width 
correction [12\ 113] , For the statistical error estimation, Delta theorem is 
used [33]. 
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Fig. 1. The raw net-charge multiplicity distributions Au+Au collisions for three 
centrality bins 0-5%, 30-40% and 60-70% and collision energy for v /snn=7.7 to 200 
GeV. 



3. Results 

The raw event-by-event net-charge multiplicity distributions are plotted 
in Fig. [T] for three centralities, 0-5%, 30-40% and 60-70% for above seven 
colliding energies. It is observed that the width of net-charge multiplic- 
ity distributions increases with increase in centrality and mean of the net- 
charge multiplicity distributions shift towards positive side with decreasing 
colliding energies. In Fig. [21 the moments such as, M, a, S and k for the 
net-charge distributions are plotted as function of (A r par t)for the above seven 
colliding energies. The M and a values increase in going from peripheral 
to central collisions for the above seven colliding energies, whereas S and k 
values decrease with increase in collision centrality. The evolution of these 
higher moments can be understood by the central limit theorem (CLT) |13j . 
which explain the (A^ part ) (proxy of volume) dependence of these moments. 

For the non-CP baseline at the above energies, the Poisson expectation 
is studied. In this approach, the positive and negative charged particles 
distributions are assumed to be independently Poisson distribution where no 
dynamical correlation among the positive and negative charge particles are 
taken into account. Hence, the net-charge distribution is taken as a Skellam 
distribution (resultant of the difference of two Poisson distributions) which 
serve as one of the baseline for this analysis. For the Skellam distribution, 
a _ n++n_ g _ n + -n_ anc } K(J 2 _ y Here n_i_ and n_ are the means 

M n+—n— ' n++ri— ~ 

of the positive and negative charged particle distribution. The results from 
the Poisson expectations are discussed and compared below. 
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Fig. 2. The various moments, mean (a), standard deviation (c), skewness (b) 
and kurtosis (d) are plotted with respect to (-/V part ) for Au+Au collisions from 
V / %n~7\7 to 200 GeV. The dotted lines represent the expectation form CLT. 
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Fig. 3. Sa are plotted as a function of (N palt ) for Au+Au collisions from v / s NN =7.7 
to 200 GeV. The dotted lines represent the Poisson expectation from the data. 



In Fig. O and HI the values of Sa and na 2 are plotted as as function of 
(-^Vpart)for the above seven energies and corresponding Poisson expectations 
are also plotted for each energy. The deviation between Sa and that of Pois- 
son expectation are observed to be increasing as decreasing colliding energy. 
Whereas in na 2 , deviation from Poisson expectations are observed almost 
similar in all energies. In Fig. [5j the energy dependence of the product of 
moments such as, fj, Sa and na 2 are plotted along with Poisson expec- 
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Fig. 4. kct 2 are plotted as a function of (N part ) for Au+Au collisions from v / s NN =7.7 
to 200 GeV. The dotted lines represent the Poisson expectation from the data. 
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Fig. 5. The jj (a), Sa (b) and no 2 (c) are plotted as a function of the colliding 
energy for three centrality bins 0-5%, 30-40% and 60-70%. The blue shaded re- 
gion represent the Poisson expectation, reflects the range covering 0-5% to 60-70% 
centrality. The red dotted lines are the HRG predictions. 
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tation and Hadron Resonance Gas (HRG) model predictions. The values 

2 

of j-j increases with increase in colliding energy for three centrality bins 
whereas in case of the Sa, it increases with decreasing colliding energies. 
The Ka 2 shows no energy dependence and all the values are above unity, 
except for most central events at 7.7 GeV with large statistical uncertainty. 
The energy dependence of the Sa shows systematically large deviation from 
its Poisson expectation below - v /s NN =27 GeV as compared at higher ener- 
gies. The HRG predictions of Ka 2 are very close to data whereas that of Sa 
over-predicts the data. 

4. Summary 

The higher moments of the net-charge multiplicity distributions have 
been measured for Au+Au collisions at \/s NN = 7.7 to 200 GeV. The cen- 
trality dependence of the moments follows the expectation from the CLT. 

2 

The values of j^- increase with increase in colliding energy. The values of 
Sa increase with decreasing colliding energies and deviates from Poisson 
expectation below \/s NN = 27 GeV. Within statistical uncertainty, Ka 2 is 
seen to be independent of collision energy. 
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